Giant positive magnetoresistance and localization in bosonic insulators 
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We study the strong localization of bosons, motivated by several experiments that suggest a 
bosonic superconductor-insulator transition in strongly disordered films. In the insulator, unlike 
for fermions, nearly all scattering paths between low-energy sites contribute positive amplitudes, 
which interfere constructively in the hopping matrix element. The localization length of bosonic 
excitations shrinks as the constructive interference is suppressed by a magnetic field. This entails 
a giant positive magnetoresistance, opposite to the analogous effect in strongly localized fermions. 
Within the forward scattering approximation, the localization length is found to increase with 
decreasing energy. Applying the same technique to hard core bosons on a Bethe lattice with large 
connectivity, we find the superfluid to emerge without concomitant closing of a mobility gap in the 
insulator. 
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Many disordered insulators display some form of vari- 
able range hopping transport [l[ , which reflects the local- 
ization of the electrons close to the Fermi level @ . When 
the hopping lengths are much larger than the distance 
between impurities, transport is very sensitive to orbital 
magnetic fields. The latter suppresses the destructive 
interference among alternative paths of electrons, and 
thus weakens the localization, leading to a strong nega- 
tive magnetoresistance (MR) This is a non- trivial 
manifestation of quantum interference in disordered insu- 
lators. In this Letter we establish the bosonic counterpart 
of the effect, which turns out to have the opposite sign. 
The MR in the hopping regime thus serve as an indicator 
for the quantum statistics of the conducting carriers. 

Recent experiments in strongly disordered, supercon- 
ducting InOz films 0, 0] have shown that a magnetic 
field not only destroys rapidly the already weak super- 
conductivity [|[ , but also induces a giant MR peak in the 
ensuing insulating state, while at largest fields the resis- 
tance often drops back close to its normal state value. 
Similar peaks in MR have been reported in amorphous 
films of TiN , Bi [l(| , and in patterned films [llj , as 
reviewed in |12| . 

The giant positive MR in the vicinity of the super- 
conducting transition is intriguing. Mechanisms such as 
shrinking wavefunctions or spin blocking effects of weakly 
interacting electrons, which may play a role in semi- 
conductors [H, hardly apply to these systems [l3j]. In- 
stead, many experimental observations in transport 
HEEL liiGil and tunneling microscopy 1^, 17 1, as well 
as theoretical model studies [l8l - Eo| suggest the crucial 
importance of remnant electron pairing in the insulator, 
despite the absence of global phase coherence 2lj. In 
the presence of strong disorder, the resulting insulator is 
expected to be a Bose glass [22| , whose low energy excita- 
tions are localized by disorder, but do not exhibit a spec- 
tral gap. The transport properties of such complex in- 
teracting systems pose interesting conceptual questions, 



which are not fully resolved vet filial |24|. 

While it is natural to expect a magnetic field to in- 
crease the resistance of a Cooper pair insulator (in con- 
tinuation of its detrimental effect on phase coherence in 
the superconductor), there is no satisfactory microscopic 
explanation of the giant effects seen in experiments to 
date. Ref. [3] attributed the MR peak to the survival of 
superconducting islands, whose properties are tuned by 
the field. A Coulomb blockade was postulated to occur 
in the islands, hampering transport of electrons, but not 
of pairs. Beloborodov et al. [26j considered an array of 
normal and superconducting grains. They showed that 
the magnetic field tunes the gaps of pair and electron 
excitations, resp., with opposite trends, which predicts a 
peak in MR. Ref. 2jJ interpreted the peak as reflecting 
the first oscillation of the charging energy of a clean ar- 
ray of Josephson islands in a magnetic field. However, 
the models 2(| 27| hardly do justice to the complexity 
of the experimental systems 0] in which no well defined 
granular structure exists, and the spectral gap for pairs 
should be washed out by strong disorder. [18[ The latter 
suggests that it is rather pair localization due to disorder, 
which induces the insulating behavior. [l3| 

In this Letter we study localization, MR and meso- 
scopic intereference effects in bosonic insulators. These 
are also relevant for recent experiments with interact- 
ing cold atoms in strong disorder potentials and artificial 
gauge fields. [IH 

We consider a simple model of bosonic or fermionic in- 
sulators: a lattice of sites, which can accommodate only 
one quantum particle, assuming strong onsite repulsion. 
This is the standard Anderson model for fermions [H, 0] ■ 
For hard core bosons the model was introduced by Ma 
and Lee [29[ who considered disordered superfluids in 
terms of preformed pairs (pscudospins) . The calculations 
below can easily be generalized to grains or islands host- 
ing many particles, as long as the charge gap on typical 
grains is much bigger than the hopping amplitude be- 
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FIG. 1: In the configuration on the right, many paths con- 
tribute to the Green's function Gfj between the low energy 
sites i,j. The locator expansion yields the matrix element in 
4'th order perturbation theory as a sum over 6 paths. The 
sign of fermion amplitudes depends on the number of occu- 
pied sites on each path (shown on the left) whereas the paths 
for bosonic low energy excitations always have positive am- 
plitudes. A magnetic field suppresses that constructive inter- 
ference, while fermions display the opposite effect. 



tween grains. A similar model was studied in Ref. [30j, 
but the importance of quantum statistics and many body 
effects on localization and interference was missed. 

We consider sites i with a random onsite energy Ej, 
uniformly distributed in E [— W,W]. They are weakly 
coupled by a tunneling amplitude tij — t between nearest 
neighbors, 



H 



E 



tij(b]bi 



b}h. 



(1) 



b\, bi are creation and annihilation operators of fermions 
or hard core bosons, resp. They satisfy the commutation 
relations [bi,bj] B = 0, {b\,bj\ B = <%[1 - 2B(1 - m)], 
where [., .]b is the commutator or the anticommutator 
for bosons (B ~ 1) or fermions (B = 0), resp. In the 
presence of a magnetic field, the hopping acquires a phase 
= te~ l ^ ij , the sum of phases </>y around a plaquette 
being proportional to the flux threading it. 

We focus on the strongly insulating regime t <C W, 
where hopping transport is expected at low tempera- 
tures. A key element characterizing disordered insulators 
is the localization length, £. For non-interacting fermions 
it is well-defined as the (log-averaged) inverse spatial de- 
cay rate of single particle wavefunction amplitudes. In 
contrast, hard core bosons are inherently interacting, re- 
quiring a generalization of this concept. In the limit 
t = single particle excitations correspond to the ad- 
dition or removal of a particle on given sites. For small 
hopping t/W <C 1, these excitations adiabatically de- 
form into dressed quasiparticle excitations, which are still 
well localized in space. In fact, one expects that all low 
energy excitations remain discrete and localized in this 
limit [13 , 24 1 . The spatial properties of such quasipar- 
ticle excitations are well captured by the retarded Green's 



function (A(t) = e lHt A(0)e- lHt ) 

G^(t-t') = -iQ(t-t')([b i (t),bl(t')}B). 



(2) 



It describes the amplitude of finding an extra particle at 
site i, and at time t after adding a particle on site 0. As 
in early studies of the Hubbard model [H[ , it is useful to 
study the equation of motion 



6(t-t')S ifi ([b (0),bl(0)}B) (3) 



This is the starting point for a locator expansion in pow- 
ers of the hopping t/W 2]. It is easy to show that 



ibi(t) = [bi{t),H]=s i b i (t) - (-1) 



Bni(t) 



(4) 



where the sum runs over the neighbors of i. We are inter- 
ested in the decay of the correlation function at large dis- 
tance. In analogy to the fermionic (single particle) study 
by Nguyen et al. , we restrict ourselves to forward scat- 
tering paths to leading order in t/W . Hence, we retain 
only the neighbors j, which are closest to 0, cf. Fig.[T] To 
leading order, we can also neglect the time dependence of 
m(t) and approximate = sign(e l ) + 0{{t/W) 2 ). 

To characterize the spatial decay of an excitation of 
given energy, it is preferable to work in frequency space, 
G £o( w ) = S-oo G i?o We^'*, and to define the boson lo- 
calization length as the (log-averaged) inverse decay rate 
of Gf (uj) with distance, 



Jim ln[|G« M/G« (w 



n-r \ 



(5) 



The frequency dependence will be discussed further be- 
low. Note that the transition to the superfluid is indi- 
cated by the divergence of £(cj = 0), where the bosons 
condense into a delocalized state forming at the chemical 
potential. 

To leading order in t the above equations furnish a 
simple recursion relation for the Green's functions with 
increasing distance. Upon iteration, this forward scat- 
tering approximation (FSA) yields Gf Q as a sum over all 
shortest paths V (of length i) between sites and i, 
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(6) 
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Setting lu = Eo we find the "wavefunction" of the quasi- 
particle excitation, which is adiabatically connected to 
the boson insertion/removal at site in the non- hopping 
limit (by extracting the residue of the corresponding pole 
in Gf ). The forward scattering and the particle rear- 
rangements on a given path are illustrated in Fig. [TJ 

For fermions, Eq. J6J reproduces the result for non- 
interacting particles [2j, y] , while hard core bosons differ 
crucially in the sign of the amplitude contributed by the 
paths. The result is conveniently rewritten as 



0,0 M 
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V={] =0,...,u=i} P 
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The difference is simple to understand. In order to ob- 
serve a particle at site i, in response to inserting a particle 
at 0, all the n v = Yfj=i n j ~ X^jX 1 ~ s g n ( £ j))/ 2 P ar " 
ticles on the path V have to move to the next negative 
energy site closer to site i, cf. Fig. [TJ Upon removing 
the particle at site i, a ring exchange of n-p particles 
has been carried out in the ground state, which causes 
the sign difference (— l) nv between bosonic and fermionic 
amplitudes. This difference has important consequences 
and will shed new light not only on strongly localized 
insulators, but also on the approach to superfluidity. 

Most importantly, Eq. ([7J shows that for low energy 
excitations, lu — > 0, in the absence of a magnetic field, 
all paths contribute a positive amplitude and therefore 
interfere constructively, unlike fermions. This difference 
manifests itself in completely opposite response to a mag- 
netic field. It is well known that hopping fermions expe- 
rience a negative MR due to the suppression of destruc- 
tive interference 0, 0). In contrast, we find the MR of 
bosons to be strongly positive, as the phases in the hop- 
ping amplitudes reduce the constructive interference of 
paths connecting the sites relevant for transport. The 
replica scaling arguments of Ref. 4] , which maps the for- 
ward scattering problem to directed polymers, should ap- 
ply also to the bosonic case. Here it predicts a negative 
perturbative correction to the hopping probability scal- 
ing as Brj^p, where ri lop is the hopping distance. At 
larger magnetic field, one finds a reduction of the boson 
localization length, whose effect on the MR is exponen- 
tially amplified in the hopping regime. 

The opposite interference behavior of bosons and 
fermions is very likely a key element for understanding 
the giant MR peak in disordered films with remnant pair- 
ing. As long as the magnetic field does not destroy the 
localized Cooper pairs, it mainly reduces the localization 
length of the pairs. Upon destruction of the pairs, e.g., by 
the Zeeman effect, the predominant carriers are fermions, 
for which a negative MR due to increased localization 
length is predicted [3, Q. Once the latter becomes large, 
the physics of loops and weak localization (neglected in 
FSA) is likely to play a role in the negative MR, as well. 
The effects of Coulomb interactions enhance the negative 
MR even further [3^ . 

Let us now discuss how bosonic insulators differ from 
fermionic ones. A beautiful series of experiments on peri- 
odically patterned Bi films has shown that Aharonov- 
Bohm oscillations in the MR start with an up- rather 
than a down-turn of resistance at low fields. This was in- 
terpreted as a signature of bosonic carriers, for which we 
provide theoretical support here. Further, the enhance- 
ment of forward scattering due to constructive interfer- 
ence implies that hard core bosons have a larger local- 
ization length than fermions, when subject to the same 
hopping and disorder. Moreover, Eq. predicts that 
has a non-trivial energy dependence around u) = 0, 
reaching a maximum at the chemical potential lj — 0, as 



is easily confirmed numerically. In fact, the presence of 
other bosons enhances the derealization tendency of an 
extra particle, whereas non-interacting fermions are in- 
sensitive to the presence of the Fermi sea. On the other 
hand, at high energies bosons tend to behave like non- 
interacting particles, as paths through occupied sites be- 
come negligible. 

We point out, that the discussion of is based 

on the FSA, which yields a recursive relation between 
Green's functions at the same ui, like in a non-interacting 
problem. However, this "elastic" behavior breaks down 
when subleading loops are taken into account and the 
interacting nature of the problem becomes manifest. In- 
sofar, the frequency dependence of £(w) makes sense a 
priori on intermediate length scales, where the FSA domi- 
nates the propagation amplitude. Instead, at large scales, 
defined by Eq. (JSJ) , must always be an increasing 
function of uj. This phenomenology is closely analogous 
to that of phonons. Their localization length increases 
with decreasing energy, which is well-defined as long as 
inelastic scattering can be neglected. 

The reduced localization tendency in low energy FSA 
is in conflict with predictions of Refs. 13|, [24j, where the 



presence of hard core bosons was argued to impede the 
propagation of an injected boson. While the arguments 
of [13j apply to a distinguishable extra particle, they ne- 
glect the exchange effects of identical particles, which in 
fact enhance the propagation. Similar ideas may have 
inspired the Bethe lattice study in 24J. Their perturba- 
tion series in t for (GfJ(w)) 2 was restricted to the subset 
of intermediate one particle excited states, yielding path 
amplitudes as a product of l/(|ej| — u>) instead of ex- 
pression Eq. ([7]). Since this subtle difference changes the 
results qualitatively, we use ([7]) here to revisit the analy- 
sis of Ref. [24| for hard core bosons on a Bethe lattice of 
large connectivity K. The latter ensures that the FSA 
is well-controlled up to very close to the superfluid. On 
a Bethe lattice, only a single shortest path connects two 
sites, interference phenomena being subleading in t, and 
statistics is irrelevant, since no particles are exchanged. 
Indeed, one finds the same localization properties for 
hard core bosons as for fermions [HI], as characterized 
by (G^(w)) 2 at large distances. 

Neglecting loop corrections to the FSA, the super- 
fluid transition occurs when E i ,dist(i,o)=ii-v o G mi u = 
0) diverges, which happens at the same (t/W) c = 
0(l/K\n(K)) as the derealization of single parti- 
cles [33j]. This is ensured by the fact that the sums 

Ei,di s tfi,o)=fi-^oo( G S( w = °)) a O = M) are both dom- 
inated by the same rare paths, which can be regarded as 

replica symmetry breaking [2|, I24JI . In the insulator we 
find results which deviate from [24j: we find that G^(cj) 
never decays more slowly than G^(0), suggesting that 
excitations at higher energy are also localized. Since at 
large K collisions of single boson excitations are very in- 
efficient in enhancing derealizations, as argued in 24J, 
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the above might imply that the insulator is fully local- 
ized even at extensive energies. We believe, however, that 
this abrupt jump from supcrfiuid to " superinsulator" is 
an artifact of the large K limit. 

This and other features are indeed expected to be mod- 
ified at finite connectivity K, and certainly in finite di- 
mensions, where the interference of multiple paths is cru- 
cial. However, one aspect of the FSA analysis at large 
K, is expected to have more general validity. The above 
suggested that the superfluid transition occurs as a de- 
localization phenomenon at oj = 0, without concomitant 
closing of a mobility gap. Since effects beyond FSA seem 
to reinforce this feature, we conjecture it to be quite gen- 
eral at disorder driven SI transitions. This contrasts with 



the criticality of the mobility edge pro pos ed in [13|, [24 1 
and similar early ideas by Hertz et aL[34j, while being 
consistent with the conjectures of Ref. |23j .We emphasize, 
however, that the Bose glass adjacent to the superfluid 
is in general not fully localized at all energies of order 
0(1). Instead, it may possess a finite, but non-critical 
mobility gap, which can give rise to activated transport 
as discussed in [l3[ . An exactly solvable model with such 
properties will be presented elsewhere [36| . 

The bosonic locator expansion may help to under- 
stand how bosons can escape localization in two dimen- 
sions, while fermions generically localize in the absence of 
spin orbit interactions. The usual argument proceeds by 
showing that a superfluid is stable to weak disorder [35| ■ 
Our approach complements this view from the insulating 
side: At low energies, all paths interfere constructively 
(which is a precursor of the establishment of a global 
phase in the superfluid). Hence, the special role played 
by time reversed paths for fermions seems negligible for 
bosons, since most paths interfere constructively anyway. 

In conclusion, we have shown that strongly disordered 
bosons exhibit localization properties and response to 
a magnetic field which is completely opposite to that 
of fermions. It would be interesting to look for signa- 
tures of the predicted localization features, as reflected 
by the dependence £(cj, B)). Apart from superconducting 
films, pro mising experimental systems are cold bosonic 
atoms [28( , where a magnetic field can be simulated by a 
rotating frame. We hope that the locator expansion may 
be pushed to approach critical phenomena and fractality 
at the superfluid-insulator transition. 

I would like to thank S.V. Syzranov, and especially V. 
Kravtsov for many stimulating discussions. 
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